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Abstract 

■ On the basis of simple principles we derive and investigate the equa- 

' tions of relativistic plasma magnetohydrodynamics (MHD) in an arbitrary 

^ , gravitational field. An exact solution describing the motion of magnetoac- 

tive plasma against the background of the metric of a plane gravitational 
wave (PGW) with an arbitrary amplitude is obtained. It is shown that in 
strong magnetic fields even a sufficiently small amplitude PGW can create 
a shock MHD wave, propagating at a subluminal velocity. Astrophysical 
O , consequences of the anomalous plasma acceleration are considered. 

bX)! 1 Introduction 

In [14] the effect of PGW on plasmalike media was investigated by the methods 
^ . of relativistic kinetic theory in the approximation when the back reaction of 

^f) ' matter on the PGW is negligible: 

\o ' 

OO ■ e < (1) 

O . 

where w is the characteristic frequency of a PGW, s is the matter energy den- 
sity {G = h — c — I). These papers have revealed a number of phenomena of 
interest, consisting in the induction of longitudinal electric oscillations in the 
plasma by PGW. In spite of the strictness of the results obtained in [1.4], the 
effects discovered in these papers have very little to do with the real problem 
of GW detection. Moreover, the above calculations show lack of any prospect 
for GW detectors based on dynamic excitation of electric oscillations by grav- 
ed ' itational radiation. There are two reasons for that: the smallness of the ratio 

(m^G/e^) = 10"'*'^ and the small relativistic factor {v^) jt? of standard plasma- 
like systems. The GW energy transformation coefficient to plasmatic oscillations 
is directly proportional to a product of these factors. 

However, the situation may change radically if strong electric or magnetic 
fields are present in the plasma. In Ref. [5], where the induction of surface 
currents at a metal- vacuum interface by a PGW was studied, it was shown that 
the values of currents thus induced can be of experimental interest. In [6], on 
the basis of relativistic kinetic equations, a set of MHD equations was obtained, 
which described the motion of collisionless magnetoactive plasma in the field of a 
PGW of an arbitrary magnitude in a drift approximation and it was shown that, 
provided the propagation of the PGW is transversal, there arises a plasma drift 
in the PGW propagation direction. The set of equations obtained in [6] is rather 
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complex and unwieldy: it is a set of nonlinear partial differential equations. In 
[7], however, it was shown that, provided the plasma is originally electroneutral 
and uniform, the solution of the above set of equations is strictly stationary, i.e., 
it depends only on retarded time. This fact permits us to substantially simplify 
the problem and to find its exact solution, possessing a number of remarkable 
peculiarities. 

2 The conditions of magnetic field embedding 
in the plasma 

As pointed out above, in [6], on the basis of a selfconsistent set of coUisionless 
kinetic equations and the Maxwell equations (i.e., on the basis of general rclati- 
vistic Vlasov equations [8]), a set of MHD equations describing the motion of 
magnetoactive plasma in the field of a PGW, was obtained. This set of equations 
is obtained in the so-called drift approximation, i.e., in the first approximation 
in the small parameter ^ : 



where lub = eH/rrieC is the Larmor frequency. However, the equations obtained 
in [6] are applicable only in the case of a strictly transverse PGW propagation, 
where the original magnetic field is perpendicular to the GW propagation di- 
rection. It is not difficult to verify that if the conditions of the strict transvesity 
of PGW propagation are not met, the equations of [6] violate the energy and 
momentum conservation laws. For our purposes it is necessary to consider a 
more general case, so in what follows we shall obtain the MHD equations on the 
basis of other principles. 

It is not difficult to sec that a consequence of the MHD equations from [6] 
is the magnetic field embedding in the plasma (MFEP) . This reflect the general 
nature of magnetoactive plasma provided that the condition (2) is met. There- 
fore, in order to describe the motion of the plasma in a drift approximation, it 
is simpler to demand at once that the MFEP condition be met. Mathematically 

this requirement means a coincidence between the timelike eigenvectors of the 

p 

plasma energy- momentum tensor (EMT), Tik, and that of the electromagnetic 
/ 

field, Tik, i-e., according to Synge [9] , a coincidence of the dynamic velocities 
of the plasma and the electromagnetic field: 




< 1 



(2) 



p 



fc 



(3) 




£hU, 



(4) 



where 



{u,u) = g, 



Hku'u'' = 1. 



(5) 
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We shall consider in this paper the EMT of the plasma as that of a perfect 
isotropic fluid 

T'^^ [e+pyv^ ^pg'^, (6) 

where 

{v,v) = l, (7) 

and e and p are the fluid energy density and pressure connected by a certain 
equation of state: 

P = p{e). (8) 

Thus is the timelike eigenvector of the plasma EMT (w' = u*), while is 
p 

the eigenvalue T**^ [sp — e), and there remain the conditions 

f 

Tlv^=eHv\ (9) 

Thus ^ are precisely the conditions of magnetic field embedding in the 
plasma (MFEP). It is our purpose to clarify all the restrictions imposed by 
these conditions on the Maxwell tensor Fik. Using the plasma velocity vector 
11*, we shall introduce the vectors of the electric field Ei and magnetic field Hi 
as observed in the frame of reference (FR) comoving with the plasma [10]: 

E.^v'^Fkf, H.^v'^Fk,, (10) 

* 

where Fki is a tensor dual to the Maxweh antisymmetric tensor Fki'- 

Fh^ \vkamF'-^, (11) 

and rjkiim is the covarianly constant discriminant tensor (see, e.g., [9]) satisfying 
the identity 



Ihklm'l — "klm ^ 



5% 51 5^^ 



(12) 



Due to (10), the vectors E and H are spacelike and orthogonal to the velocity 
vector: 

{v,E)^0- iv,H)^0. (13) 

The relations (10) can be resolved with respect to the Maxwell tensor [10]: 

Fik = ViEk -VkEi- Tjikhnv' H'"^ ; 

Ak^ V^Hk - VkH, + Tl.kUnv'E"'. (14) 

Let us represent the EMT of the electromagnetic field as follows: 

^{F\F\ + i<5^F'"F;„0 (15) 



3 



using three vectors {v;E;H), one of which, (w), is timeUke, while two others, 
{E and H), are spacehke: 
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k- 
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+ 2E'Ek + 2WHk + 2v'TjkpqsEPv'iH' + 2v''t]'P''' EpV^Hs] , (16) 
where the following notations are introduced: 

E^ -{E,Ey, -{H,H). (17) 

We shall require that the vector v be an eigenvector of the EMT (16). Con- 
tracting (16) with and taking into account the identities and p^ . we 
get: 

^[v\E^ + H^) - 27fP'^'EpV,Hs] = SHvK (18) 

From ((T8)) we shall find the necessary condition for being an eigenvector of 
the electromagnetic field EMT: 

rj'P'i'EpVqHs = \v\ VA G R. (19) 

Contracting this relation with u*, we get A = 0. Thus the necessary condition 
of compatibility ^ and (g]) is: 

r]'P^'EpVgHs - 0. (20) 



The necessary and sufficient condition for the fulfilment of Eq. ((20| is, as we 
know, the complanarity of the vectors E, H, v, i.e., 

av^ + I3E^ + -fH, = 0. 

Contracting this relation with and taking into account (O, (fT3|) we obtain 
a = 0. Thus, Eq. ([^0]) is equivalent to the condition: 

^E, + -iH, = 0. 

Since we consider magnetoactive plasma, we shall further assume: 

F^kF'^^ = 2{H^ - E^) > 0. (21) 
Due to ([21]) . the necessary and sufficient condition for the fulfilment of ([20l) is: 

E, = XHf, VA e R (22) 

Besides, according to ([TOl) , 
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However, we have not yet extracted all the algebraic information contained in 
([3]) and (|4]). With due account of the definitions of the vectors E and H (10), 
Eqs. pU)) can be written as: 

T.kv'' = O, (24) 
where a new antisymmetric tensor J-'ik has been introduced: 

J'ik = Fik — A Fik ■ (25) 

The relations (1^^ can be regarded as a set of linear homogeneous algebraic 
equations with respect to the velocity vector v^. The necessary and sufficient 
condition for a nontrivial compatibility of these equations is: 

Det||J-,fe|| =0. (26) 

As is an even order antisymmetric matrix, 

1 ' ■ — ifc * \2 



Detll J-,fe|| = —i^r" J^^k). (27) 
Therefore, the condition p6|) reduces to the following: 

V^'^'F^jJ^ki = 0, (28) 

- in this case 

rank||J-,fc|| =2, (29) 
i.e., the set (|24p admits two linearly independent solutions for the eigenvector 

Substituting into ([28]) the expressions for Fik and Fik from (fT4l) . we get: 
v'^'^'F^.Th = 4A(1 + A2)ff2 ^ 0, 
hence follows the only possible solution under the condition ([21]) : 

A = 0. 

Thus we come to the following rigorous conclusion. For the EMT of electro- 
magnetic field (15) to permit as its eigenvector the dynamical velocity vector of 
the isotropic perfect fluid under the condition (21 ), it is necessary and sufRcient 
that the electric field intensity vector in the comoving FR should be equal to 
zero: 

E, - 0. (30) 
In this case the conditions ([24)) - ([29|) give: 

i^'" F,k - 0; (31) 
Det||i^,fc|| = =^ rank||F,fe|| = 2, (32) 
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while the eigenvector of the fluid EMT must satisfy the set of linear homogeneous 
algebraic equations: 

F.kv'' = 0. (33) 

Note that if Eq. (15^ is fulfilled, then similar conditions for the dual Maxwell 
tensor are automatically fulfilled as well: 

Dctll F,k 11=0=^ rankll F,k \\ - 2. (34) 
With (|5D|) write down the Maxwell tensor and the electromagnetic field EMT 



Fik= ViHk -VkHf, (35) 

^ f 

rfc= —i2H\'vk - 2WHk - 51H^). (36) 
an 



Note that due to Eq. ([35|) a more severe condition than (|3ip is fulfilled: 

F,k F"' = 0. (37) 
The summed EMT of the magnetoactive plasma 

_p f 
Tik ^Tik + Tik 

takes the form: 

T^k = {£ + P)viVk - Pgik - 2PHnink, (38) 

where 

Ph = ^; 8 = e + eH] P^p + Ph, (39) 

OTT 

P and £ being the summed pressure and energy density of the magnetoactive 
plasma and 

n. = § (40) 
- is the spacelike unit vector of magnetic field direction: 

{n,n) = -l, (41) 

with 

(n, v) = 0. (42) 
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3 MHD equations for plasma in a gravitational 
field 



The MHD equations are to be obtained on the basis of the vanishing divergence 
requirement for the summed EMT of the magnetoactive plasma (j38p supple- 
mented by the first group of the Maxwell equations: 

= 0, (43) 

F'\k^ 0. (44) 

This set of equations with due account of the equation of state (8) , the definition 
of the plasma EMT (0) and the algebraic relations 0, (HH), dSO]), ([35]) 

and ([55)1 completely describes the self-consistent motion of the magnetoactive 
plasma with an embedded magnetic field in a prescribed gravitational field. 
Indeed, Eqs. dU]) represent a set of 8 differential equations with respect 

to 10 quantities £,p, However, the equation of state the velocity 

vector normalization (O, the ortogonality condition (IT3l) and ([30]) raise the 
total number of equations up to 12. Nevertheless, it turns out that not all of 
these relations are independent, as we shall see below. 

To deduce the MHD equations, let us take into account the well-known 
relationship (see, e.g., [11]): 

T'^,k^ -^F^iF'i- (45) 
Thus, Eq. can be presented in the form: 

F.fc*'^ = T„ (46) 

where 

r,; -An T^,fc . (48) 

Eqs. pSj) can be regarded as a set of linear inhomogeneous algebraic equations 
with respect to If Det||i^ifc|| ^ 0, then the equations are solved in a simple 
straightforward way: 

(49) 

Det||i^|| ' ^ ^ 

where A''\\F;t\\ is a cofactor of the augmented matrix of the set (|46p . In partic- 
ular, for the case of vacuum (t^ — 0) we obtain a trivial solution: i"*^ = 0, and 
the set of equations (|^ . (El) reduces to the Maxwell equations in vacuum. 

Now we pose the problem of solving the set of equations (HH) with respect 
to i.e., the problem of reducing Eqs. (l43l (|44)) to the form of the Maxwell 
equations with minimal requirements upon the electomagnetic Geld invariants: 

F''' hk= 0, (50) 
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F^kF''' > 0. (51) 

The existence of a positive invariant (|5ip means that we can choose a local FR 
where the electric field is absent [11]. 

Provided (50) is fulfilled, as before (see (|32|) ). 

Det||F||=0=^ (52) 

rank||i^|| = 2. (53) 



Thus for the consistency of the algebraic set of equations p6|) under the condi- 
tion (|50[) it is necessary and sufficient that: 

rank||F;T|| = 2. (54) 

Calculating all the 3rd order minors of the augmented matrix r|| with (jSO]) . 
we obtain a condition equivalent to 

hk Tk = 0^ rank||i^; r|| = 2. (55) 



To solve the set of equations (|46l) . consider the eigenvectors of the matrix 

mk\\: 

F.ku'' = Xu,. (56) 

Due to the antisymmetry of Fik, it follows from that either A = 0, or u 
is a null vector. It can be demonstrated that provided Eqs. (|50p and ([5T|) are 
valid, M cannot be a null vector. Thus, the Maxwell tensor admits only nonnull 
eigenvectors with zero eigenvalues: 

F^ku'' = 0. (57) 

By ([50)l . ((52|) and ([SS]) . this eigenvalue is doubly degenerate, and thus, according 
to a well-known algebraic theorem, two lineary independent eigenvectors cor- 
respond to it: u and u . 

(1) (2) 

Under the condition ([5T|) we can always choose a local FR where Fa — 0. In 
this FR there always exists an eigenvector of the Maxwell tensor of the form = 
5^. Therefore, one of the eigenvectors of the matrix ||i^||, e.g., u , is timelike and 

the second one, u , is spacelike. Using the standard orthogonalization process, 

(2) 

we normalize them as follows: 

(u,u) = l; (u,u) = -l; ( u , u ) = 0. (58) 

(1) (1) (2) (2)' \l) (2)' 

Then the general solution of (j57l) can be written in the form 

^au^ ^I3u^, (59) 
(1) (2) 

where a and /3 are arbitary scalars. 
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Let us now investigate the relations (I55|) , which can be regarded as algebraic 
equations with respect to r. Since 

- F.feF^'^ ^hkF'''< 0, (60) 

and Eq. ([50)) is invariant under the substitution F ^F, as well as the expres- 

* 

sion for Det||F|| ([?f|) . we conclude that the dual matrix || Fik \\ also admits 

two and only two linearly independent spacelike eigenvectors w and w , which 

(1) (2) 

correspond to a zero eigenvalue. It is not difficult to verify (e.g., turning to a FR 

where Fai — 0) that the rank of the unified matrix \\F, F \\ under the condition 

(|50|) is equal to 4. Hence the eigenvectors of the matrices ||i^|| and || F || are 
linearly independent and we can choose the following normalization for them: 

(1) (1) (2) (2)' 

(«; u;)=0; C^,.") = 0> (a,/? = 1,2). (61) 

(1) (2) (a) (/3) 

Thus, the general solution to Eq. ([55)) is: 

Ti = A +^ w^, (62) 

(1) (2) 

where A and n are arbitary scalars and due to ([61)): 

{tu) = 0. (63) 

(a) 

By ([5^ and (1551) . the Maxwell tensor and its dual can be represented in 
terms of the eigenvectors of the matrices ||F|| and || F ||: 

Fik = -uVikim u^u""; (64) 

(1) (2) 

F,k^ Q-qiklm (^^'^^j™' (6^) 

where a and g are certain scalars. Contracting these relations with the discrim- 
inant tensor, we obtain the dual relations: 

F'^=a(u'u'' - u'u^)- (66) 

\l) (2) (2) (1) ' 

F^'' -piw'w^ - w^w^). (67) 

(1) (2) (2) (1) ' 

In particular, the following relation stems from ()64p and (|66p: 

Ak F^' = 0. (68) 
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Using the Maxwell tensor representation ([M)) and the orthonormality rela- 
tions ([58|. we get the useful formula: 

F^kF'' = a^i- Uku' + Uku' +4). (69) 

(1)(1) (2)(2) 

Contracting we get: 

^F.kF''' = a2 > 0. (70) 

Contracting Eqs with F'' and taking into account and (1701) , we obtain 
an equation equivalent to Eq. (j46p : 

cr2[$'_ (w ,$)+ (m ,<I>)1 = F'^i. (71) 

(1) \l) (2) \2) 

A special solution to Eq. ([7T|) is: 

= (72) 

Therefore the general solution of Eqs. can be presented in the form: 

p 

= = ^ ^ + a +/3 (73) 

■'^ F'™F„ (1) (2) ^ ^ 

This exhausts the problem of reducing the set of equations (H5|) . p^ to the 
standard Maxwell form. 

If the external currents are absent {a = (3 — 0), then Eq. ([75)) reduces to 
the form of the second group Maxwell equations: 

= -Ml, (74) 

where: 

p 

is the drift current^ which by (j68p satisfies the relation: 



Hence, due to ([55)1 and 



(•^rfr'A)=0, (77) 

(q) 

(Jdr,Jdr)<0- (78) 

i.e., i/ie drz/t current is spacelike. 

Calculating the covariant divergence V/ of Eq. with the aid of the first 
group Maxwell equations (|44|. we get the differential implication: 

F''' Fku = F""' Fku= 0. (79) 
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Calculating the covariant divergence V of using the second group Maxwell 
equations (|74p and the relation (j76p . we get one more differential implication: 



F^' F,ik,i]=0. (80) 

Setting in Eqs. 1^,1^, §^,§^-§7!) v\ u'^ = (H'/H), we 

find a complete coincidence of the above formulas with the relevant expressions 
from the previous section. Thus: 

= = -{H, H) = H\ (81) 

Thus due to (|63| the following differential relations should be valid: 

v' f\^k= 0, (82) 

H' T^,fe= 0. (83) 

(pl the 

(O and (US]), we obtain: 



Substituting to ([82)) and (|83)) the expression for Tik ^ and taking into account 



v% - (84) 

v,,kH^v'' = (85) 
Calculating the drift current ((TS]) with Eqs. ©,([33]) and ([51])), we get: 

J^ ^ F'''[vk,iv\s+p)-p,k] . . 

-^dr - p pim ■ W 

From f[5F)) we obtain the useful relation (for H ^ 0): 

= (87) 

By (35) and the definition of the vector W (10) the orthonormality relations 
([13]) are fulfilled identically. 

Note that due to ([68|) the solution to Eq. ([35]) is 



=F''' Sk, (88) 
where Sk is an arbitary spacelike vector, satisfying the only condition: 



F'^'F^ SuSi = 1. (89) 
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Now let us turn to the first group of the Maxwell equations (|44l) . Using the 
representation ([35]) for the Maxwell dual tensor, we obtain for Eq. (|44)) : 



v'H^^ + v'j^H'' - v%W - v^H'j^ = 0. (90) 

Consecutively contracting Eqs. (|90p with and i?' and using ([7|) and p^ . we 
obtain: 

- w.^fciJ't;"' = i?.,fei^'w'= - (91) 

H^^kv'H'' = -v^.kH'H'^ = H{Hv''),k- (92) 

* 

Since the rank of the matrix || Fik \\ equals 2, the relations (1^0]) - are 

equivalent to the first group of the Maxwell equations (H^ . By (|?T|). Eq. (|85p 
reduces to a form similar to ((M)) : 

4 Solution of magnetohydrodynamic equations 
against the background of a plane gravita- 
tional wave 

4.1 Initial conditions and symmetry of the problem 

The metric of a PGW with the polarization e+ is described by the expression 
[12]: 

ds^ = 2dudv - L^[e^^{dxY + e-2/5(dx3)2], (94) 

where /3{u) is an arbitrary function (the PGW amplitude); the function L{u) 
(the PGW background factor) obeys an ordinary second-order differential equa- 
tion; u — -^{t ~ is the retarted time and v = -^{t + x^) is the advanced 
time. The absolute future is represented by the region r+ : {u > 0;v > 0}, 
the absolute past by T~ : {u < 0;v < 0}. The metric ([M)) admits the group of 
motions G^, associated with three linearly independent Killing vectors: 

e^Sl; e=Sl; e=Sl (95) 

(1) (2) (3) 

Let there be no GW at w < 0, i.e., 

^(")|«<o = 0; L(w)|„<o = 1, (96) 



the plasma is homogeneous and at rest: 

1 

£\u<Q = ^o; P\u<o = Po, (97) 
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and a homogeneous magnetic field vector belongs to the plane {x^, x'^}: 
Hi\u<o ^ Hacosfl; H2\u<o = Hosinfl; 

H3\u<o = 0; £^j|u<o = 0, (98) 

where fl is the angle between the axis Oa;^ (the PGW propagation direction) 
and the direction of the magnetic field H. The conditions agree with the 
vector potential: 

Ay = A„ = ^2 = 0; 

A^^ Hoix^smfl-x'^cosn); {u<0). (99) 

In [7] it is demonstrated that the solution to the MHD equations in the metric 
([M)) with the initial conditions ([M]) - (IM)) with e ^ is strictly stationary, i.e., 
all observed quantities are functions of solely the retarded time u. Therefore, we 
shall immediately require that the solution of our problem inherit the symmetry 
of the metric ([MD: 

£ P = 0; (a = 1,2, 3) (100) 

for all observed quantities P (£ is a Lie derivative in the direction ^): 

p^p{u); £ = £(u); = w*(m); (101) 

F,k = F,k{u). (102) 

With (|102p we obtain the following results from the first group of the Maxwell 
equations: 

^2 pua^ ^^^^ Const); a = {w, 2, 3}. (103) 

Thus, using the initial conditions ([98]) - (|99l) . we find: 



^2 puv^ _p^_^ ^ HncosVt; 



The condition ^ {F,k F"-^ = 0) with (fT04l) reduces to: 

j2 pv-i^ p^^ ^ \piFuv cot f7, (105) 
while the MFEP conditions with pIM)) and yield: 



3 V^Fuv 



Hq smi2 

-^HoVu sinfi + FusVy — Hqv^ cos^l — 0. (107) 
v2 
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4.2 Conservation laws and the vector potential 

A consequence of the second group of the Maxwell equations (fM]) is, as we know, 
the current conservation law, which, in view of (jlOOp and the initial conditions 
dnZI) - dMl), takes the form: 

Jdr = 0. (108) 
Calculating J^,, using §6^, (pM]) and (IT06| . we reduce ([T08)) to the form: 

v^iL^e-^^F^J = 0, (109) 

where a prime is a derivative with respect to u. Thus, due to the initial con- 
ditions {Fuv{u)\u<o — 0) and the the velocity vector timelikeness requirement 
{vv ^ 0) we draw the conclusion that the current conservation law (|108p is 
equivalent to the requirement: 

Fuv = 0. (110) 

But then by piI5| and pUg]) 

F'^'= Fu2 = 0; (111) 

= 0. (112) 

As we know, the first group of the Maxwell equations (|44p is equivalent to the 
existence condition of a vector potential Ai: 

F,k^d^Ak-dkU (113) 

Thus, we can write for the zero component of the Maxwell tensor (|104p . (jllOp . 
(EID: 

d^A^-d^A^^O- ({7,a} = {m,i;,2}). (114) 

As known, a unique solution to Eqs. (jll4p on the 3-dimensional hypersurface 
E3 : {x^ ~ Const} is a gradient vector: 

A, = {a = u,v,2), (115) 

where $ = $(u, x^, x"^) is an arbitrary scalar function. 

The nonzero components of the Maxwell tensor, F^^, can be represented, by 
ifTTS]) . as: 

i^a3 = a,i3, {a^{v,u,2}) (116) 



where 



is A3 - 93$ (117) 



is a gradient-renormalized vector potential. Calculating the nonzero components 
of the Maxwell tensor F^rs using (|116p - (|117p and taking into account the 
relations pH^ and piH)) . we finally find: 

^3 = -Hox^ cosn + -^Ho[v - ^{u)]smVL, (118) 
v2 



14 



where ip{u) is an arbitrary differentiable function satisfying the initial condition 



V'|u<o = '"- (119) 

It should be noted that the only '(/'-dependent nonvanishing component of the 
Maxwell tensor is 

F„3 = -i=iJo^'sinl], (120) 
and the MFEP condition pIT7| takes the form: 

^(v^V'' - w«) sinl7 + cosfi = 0. (121) 
v2 

Calculating the other components of the drift current using ( (|104p . (|110p and 
PTT|) . we find: 

r = J^ = 0, (122) 

and for the only nontrivial component we obtain an expression coinciding 
with the one found in [6] in the case of rigorously transverse PGW propagation 
(cosfi = 0). However, for other values of f2 the expression for the drift current 
obtained in [6], as well as that for the drift velocity (cf. (|121|) ). is erroneous. 
However, we shall not integrate the Maxwell equations with the drift current, 
since it is a consequence of the conservation laws (|43| and the MFEP conditions: 
it is much simpler in our case to integrate the conservation laws directly. 

4.3 Integrals of motion 

Due to the conservation laws for the complete plasma EMT and the pres- 
ence of three linearly independent Killing vectors (|95l) there are 3 integrals of 
Eqs. dSJ: 

i'^'Tfe^C; (a -1,2,3), (123) 

(q) (a) 

where is described by Eqs. §^ and 

Calculating the magnetic field vector W and the scalar H2 according to (jlOp 
and ((8T|) . we find: 

Hy = —HoL~^{vy cos n ^ — pt;2sinr2); 

Hu = HoL^^{vu cos il -j=%p'v2 sin f2); 

V2 

H2 = -^HQe^f^smn{vu + VyiP'); H3 = 0; (125) 

v2 
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cos^ Q. sin Q. 



L4 



-#6 



(126) 



It is not difBcult to verify that according to (|125p the orthogonahty relation (I13p 
((f , i?) = 0) is satisfied identicahy. The velocity vector normalization condition 
(O with (I12ip and (|125l) can be written in the form: 



1 ■ "\2 _ -^\.2r4 



(w^cosilH — ^t;2sinr2) — —^v^L 
v2 -Wq 



sin^ 



(127) 



One of the integrals (1123p . C, proves to be trivial, while the remaining two yield 

(3) 

in view of dm]), ((nSl) and ^2B), dEi, (IT27| : 



(e+p)L2z;2 + 



sin^ rj 



,2/3 



(1) 



1 



^2 sin^ n 



4tt 



(e +p)L'^VvV2 



Thus: 



i?o cos sin 51 2^ 
4V27r ^ 

£0 +Po 



(2) 



ffp cos fl sin 



4V27r 



A(«); 



where: 



— = \/2(A"^ - l)cotr2, 



A{u) l-a^{e^^-l) 
and a dimensionless parameter a has been introduced: 

2 sin^ 



47r(£o +Po) 



(128) 
(129) 

(130) 
(131) 

(132) 
(133) 



Let us now turn to the consequences of the Maxwell equations. Integrating 
, we find one more integral: 



\/2L'^\vv\ = exp 



de 



-p{e) 



(134) 



Thus, if an equation of state is specified, p = p{e), dH), then using (|130l) . (I13ip 
and (|134p . the functions Vv{u),V2{u), e{u) and p{u) are determined. However, 
to be able to determine Vu{u), it is necessary to find the function tlj{u), for which 
yet one more integral is required. It is precisely the normalization relationship 
in the form (|127p that is the integral in question. Multiplying (|128p by cos fl 
and (|129p by sinfi and adding up the values thus obtained, we gell^: 



L'^Vtj{vv cos r2 



V2 



V2 sin VL) 



£0 +Po 
2(£+p) 



cos 57. 



(135) 



'It is easily demonstrated that this integral is directly obtainable from l|93|l . 
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Squaring both parts of (jl35p and taking into account the normahzation condi- 
tion in the form (|127p afong with (|130[) . we find: 

2 Hq / cos'^ fl e + p 



H' = -ir( ^TTlF + -^^e^^ sin2 n . (136) 



Comparing (I136P and (|126p . we finahy obtain: 



L'^{s+p) ^ f 1 _2^cot^f^ 



A(£o+Po) \A^ J 
Finally, using (|12ip . (|13ip and (|137p . we find an expression for v„: 

I^(£±rt ^/l ^y^_.,co^_ (138) 
vv A(£o+Po) \A J 

Thus, with a specifeied equation of state ([5]), all unknown functions are 
found, the solutions to Eqs. ( p^ . p^ . pM)) - (in5)) automatically satisfy- 
ing the initial conditions ((97)) . (|98p (see also (jll9l) V Thereby we have found 
by quadratures an exact solution to the self- consistent problem of motion of a 
magnetoactive plasma against the background of a PGW. 



5 A study of the solution 

Squaring both parts of Eq. (|134p and using, in the lefthand side of the equation 
thus obtained, the expression for from p30p . we represent (jl34p after certain 
obvious transformations in the form: 



A{u) ^ e--'^'-'°\ (139) 



where A =^ L'^{u)A{u), 



J(e;eo)- / —^de. (140) 
J e+p[e) 

With a specified equation of state ([8]) , Eq. (|139p completely determines the func- 
tions e{u) and p{u) and thus explicitly determines the solution to the problem 
posed. Let us investigate these functions, making the most general assumption 
on the equation of state: 

pie) < e. (141) 

Then 

P'e < 1, (142) 

e + p{e) 
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and thus 

sgn[ J (e; eo)] = sgn(e - £o) (143) 

It can be seen from (|139p that e depends on the retarded time only through the 
function A(m): 

e = e{A{u)). (144) 
5.1 Singularity investigation 

Let us ivestigate the dependence e(A). It fohows from (I130p that the solution is 
specified in the interval A e [0, +cx)). As A ^ +0, according to (|139p J{e; Eq) 
+00, which by (jl43l) is possible only with s — >■ +oo. Thus we can draw a general 
conclusion that the solution of the MHD equations against the background of a 
PGW contains a physical singularity on the hypersurfaces u ^ u^: 

A(u,) = 0. (145) 

By definition of the functions A(it) and A(w), there can be two types of such 
hypersurfaces: 

B). l-a2(e2'3(")-l) = 0. ^^^^^ 

The first type of singularity is well-known: it is connected with a coordinate 
singularity of the metric (|94l) and always arises in plasma (see, for instance, [1]). 
The second type is new and is not connected with a coordinate singularity of 
the PGW metric (|94|) - this is a purely physical singularity [13]. By ()146|) . the 
conditions for the formation of a second-type singularity are: 

P{u) > 0; (147) 

> 1. (148) 

It is well-known (see, e.g., [12]) that the values of /3(u) > correspond to a 
compression phase of the geodesic tube along the 0x2 axis, while Ox^, f3 < 
correspond to an expansion phase. But the condition (|148p . according to 
(|133|) . means that in the initial state Eh sin^ $7 > ^(eo + Po); i-e., the plasma is 
highly magnetized. It is an extremely important fact that the B - type singular 
state is possible even in a weak PGW (|/3| <^ 1) provided the plasma is highly 
magnetized (a^ 1); in this case, according to (|146p , the singular state occurs 
on the hypersurfaces u — u^: 



f3M = (149) 



Further, according to (jl39l) and (|140p . at A(u) = 1: e = eo,p = po and thus by 
(|130p - (|138p the initial conditions are restored on the hypersurfaces A{u) = 1. 
As A(u) +O0, by (|139l) J{e; Eq) — oo, which is possible by (|143p only when 

E +0;p{£) +0: 

A{u) +00, E^O; p{e) 0. (150) 
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Let us investigate the general behaviour of the solution in the vicinity of a B- 
type singularity, i.e., as A(m) — >■ 0, imposing a more severe requirement than 
(|14ip for the equation of state: 

p{e) < e/3. (151) 
Then the following inequalities are valid: 

2e - e + p{e) e ^ ' 

the equality in the left-hand side (I152p being achieved only in the case of the 
ultrarelativistic equation of state (s = 3p). Then, restrictions upon J(e;eo) 
follow from (fT40 l) . (fTsS ]) : 



hence by (fT39| : El 



In. /— < J(£;£o) < In— , (153) 
£o So 



lei 

A £o A2 



But then we obtain directly from (jl30|) : 

Il'A'<vI<Ia'. (155) 

Using (|154l) . (|155l) and singling out the dominant parts of the expressions 
(|13ip and (|136l) - (|138p near a B-type singularity, we obtain the asymptotic 
estimates B: 

s^eoA-^; (^^£(1,2]) 

A->0 ff^ifoA-^; (m €(!,§]) (156) 

w„K-A-^; (7e(2,3]) 
W2/wd ~ \/2 cot riA^-'- 

The components of the 3- vector of physical velocity!^" = dx'^/dt can be pre- 
sented in the form; 

yl _ Vu/Vy - 1 . 

Vu/Vv + 1 ' 

L2(1 -I- Vu/Vy) Vv 

Using the estimates (|156p for the functions near the singularity, we obtain for 
the physical velocity vector components: 

Fi^l; 

2 



0: , \,-2,i , , . (158) 

-^^cotr^A-'-i ^0 ^ ' 



^li, instead of (151), we restrict ourselves to the weaker condition (141), then only a lower 
bound upon the energy density will be preserved. 
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Thus a physical singularity occurs on the hypersurfacc u = the energy 
densities of the plasma and the magnetic field tend to infinity, the dynamic speed 
of the plasma as a whole tends to the speed of light in the PGW propagation 
direction. Meanwhile the plasma transverse motion vanishes {V'^ — 0), whereas 
the ratio of the magnetic field energy density and that of the plasma tend to 
infinity: 

A^O; — ^ a'e^^A-i ^ oo. (159) 

e 

The presence of a S-type singularity naturally poses the question of the 
applicability of the MHD plasma model in the vicinity of the hypersurfacc (|146p , 
i.e., the question of fulfilling the drift approximation condition ([5]). 

Note that ([2])) is a local condition, i.e., for the MHD plasma model to be 
applicable, this condition must be fulfilled in the comoving FR throughout the 
whole domain under consideration. Let us assume that ^ is fulfilled in the 
initial state u u < 0: 

Co - ^ « 1, (160) 

where lob— ^b\(h=Ho)- The GW frequency in a frame of reference moving at a 
speed w*, is uj' ~ kiV^, where k is the GW wave vector: k — (— a;,0, 0,w). Thus 
oj' — \/2vvUJ. The scalar H ([T7|) is the modulus of the magnetic field intensity 
in the comoving FR. So the local value of the drift parameter ^, measured in 
the comoving FR, is related to the initial value of (|160p by 

e = y2^eo. (161) 

Using the asymptotic estimates (|156p and (|158p of the solution behaviour in the 
vicinity of a B-type singularity, we obtain from (jl6ip : 

A ^ : ^ - ^qA^" 0. 

Thus we can conclude: if the drift approximation applicability condition was 
initially fulfilled, then in approaching to a _B-type singularity the plasma motion 
is more and more precisely described by the MHD model. Therefore, within the 
scope of the problem in hand, viz., that of motion of initially homogeneous 
magnetoactive plasma against the background of the PGW metric there are no 
mechanism available for preventing the singularity. 

We should note an essentially nonlinear nature of the phenomenon detected: 
if the initial equations are expanded in a Taylor series in the small PGW am- 
plitude, f3, this phenomenon is not observed, as confirmed, in particular, by 
an investigation carried out in Ref. [14]. The reason for this difference of the 
results lies in the fact that the governing function of the process under study is, 
as can be seen from the solutions to (|130p - (|132p . the function A~^, which in 
the case of a weak PGW (|/3| ^ 1) takes the form: 
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The Taylor expansion of this function in powers of j3 assumes the smallness of 
the value 2a^/3; however, the parameter a in a highly magnetized plasma may 
prove to be so great that the condition 2a^|/3| > 1 is fulfilled. 

Thus, the discussed phenomenon is a threshold effect arising when the PGW 
amplitude in the compression phase along the axis Qx^ reaches the value (|149p 
and can be interpreted as a non-linear threshold GW-generation of a magneto- 
hydro dynamic shock wave propagating at a suhluminal speed along the PGW 
propagation direction. We shall further call this new class of effects gravimag- 
netic shock waves, or GMSW ior short. 

Let us note that, firstly, the presence of a plasma qualitatively changes the 
electromagnetic field nature [7], [14]: in the action of a PGW on a vacuum mag- 
netic field, the solution of the Maxwell equations is essentially nonstationary (an 
exact solution is given in [6]) and, secondly, neither in a fluid, with any equation 
of state (exact solutions being given in [15], [16]), nor in a vacuum magnetic field 
singularities other than an A-iype coordinate singularity occur. The only and, 
at the same time, exotic example is the plasma with a scalar interaction in the 
case of repulsion of two identically charged particles (an exact solution is given 
in [17]). It should be also noted that in the case of the ultrastiff equation of 
state (p = e) the solution to the hydrodynamic equations is nonstationary, i.e., 
depends on the variables u and v [16]. The same type of behaviour is detected 
for a magnetoactive plasma: as seen from (I139P - (I140p . at p = e there is no 
stationary solution. 

The shock wave formation mechanism seems to consist in the following. A 
weak GW is known not to interact with a fluid but to perturb a magnetic 
field. This in turn causes a drift of the plasma. Particles with smaller values 
of the coordinate have a greater value of the coordinate u, that is why at 
the compression stage of the geodesic tube in the direction Ox^ : (/3 > 0; /3' > 
0) such particles have a greater speed than those with larger values of the 
coordinate x^. Therefore, the backward plasmic layers overtake forward ones 
and thereby contribute to the shock wave formation. 

Calculating the drift current with the aid of the Maxwell equations (17^ . we 
find (see also [6], [14]): 

J^-^iLFj ^/3. (163) 

For this reason, in a singular i3-type state the drift current density remains 
finite. 

5.2 Barotropic equation of state 

We shall examine a barotropic equation of state: 

p = ke,0<k<l. (164) 

Then 

(1 + fc) eo 
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and thus, according to (|140p . (|130p - (|133p . we obtain an exact solution: 

£ = eaK''^+^] (166) 

«„ = -^i'=Ai+t; (167) 

^ = [A-« + (A - l)2L-2e-2^ cot^ f]l ; (168) 

H^^^ = (cos2 r! + L^A-e^^* sin^ n) , (169) 

>0. (170) 



where 



1 - k 

In particular, for a nonrelativistic plasm£[f| (p = 0, k = 0) (|166p - (ll69p 
imply: (fT66| - (fT69l) : 

A-l 1 A 

^ - A-2 [1 + (A - l)2L-2e-2/s cot^ f]] ; 

= ^^{cos^n + L^sm^ne^l'); (171) 

and for an ultrarelativistic plasma {p ^ e/3 ^ k — 1) we obtain from (|162p - 
(UMl): 

e = eoA-^ w^ = -i=LA3/2; 

^ = A-2 U-i + (A - l)2e-2^ cot2 17] ; 

H^ = ^{cos^n + A-'e^^sm^n). (172) 

The exact solutions ()166p - (|169p confirm the asymptotic estimates for the 
general behaviour of plasma near the sin gularity ([T56l) . 

5.3 Drift of a nonrelativistic plasma in a weak PGW 

Let us examine the practically important case of a weak PGW: 

|/3(u)| « 1 (173) 

and a nonrelativistic plasma (fc = 0). It is known (see, for instance, [12]) that 
^ O{0^). In view of this fact, we obtain from ()17ip in the first approximation 
with respect to /3 (but not with respect to a^/3!): 

e = eoA~^; i;^, = -^A; 



^The exact solution for a nonrelativistic plasma was obtained in [18] 
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^ = A-2[l + (2a2/3cotf])2]; 

Vy 

H^^; — = V2{A'^ - 1) cot (174) 

A Vy 

where it is necessary to use the expression ()162p for A~^. Then we find from 
P57p the nonzero physical speed components in an explicit form: 

T/i - o„,2^ l + ag^cos2» _ 



^ -7!""^--^ ^-.^/3V"2lt/?2cos^^) ' (1^^) 

with = H^/4TT{eo+po). 

It follows from (1175^ that in the case of a sufficiently weak PGW (2a^|/3| < 
1)V^ > for /3 > and < for /3 < 0, i.e., in the compression phase of of 
the geodesic tube in the direction of the Ox^ axis, the plasma drifts as a single 
whole in the the PGW propagation direction, whereas in the expansion phase 
it drifts in the opposite direction. 

Let /3(u) be a quasiperiodical function with the period T , so that 

(/3H) = 0, (176) 
where (■ ■ • ) denotes averaging over the interval Au — u — uq ^ T ; 

|g2,3(«+T)_g2/3(n)|^^ 1 . 

i.e., the PGW amplitude is little changed in the course of the period T. Setting 
in (I175P = tt/2 and averaging V-^ under the conditions (|176p - (|177p over a 
sufficiently large interval of retarded time, we find: 

{V')^2atW'iu))^- (178) 

thus the average plasma drift proceeds in the PGW propagation direction and 
even in a weak GW it is nonzero. At V. — 77/2 by (I174p Vu/vy = — A~^, 

therefore: 

dv 1 1 



du A? (1 - 2all3{u)Y 



^"""^/(l-2ag/3(u0)- ^'''^ 



Passing on in (|179p to the u ^ u^, under the condition that the function /3(u) 
be continuous, we get: 

u~,u.; ,iu)^vo^j^^^^^j^^^oo (180) 
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Eq. (|179p described in an implicit form the trajectory x = — f{t) of an 
infinitely small volume of plasma. Let this volume be placed at a point xo 
prior to PGW arrival, then — a;o is the moment of PGW arrival at this 
point. The singularity occurs aX u = u^, = -^{t — x)^. Then it follows from 

(|180p that in this state x* +00 ;t* — >■ +00. Thus, by an external observer's 
time, a plasmic particle reaches the singular hypersurface u = u^, (which the 
plasma velocity tends to that of light and the energy density tends to infinity) 
for an infinitely long time, and, in doing so, it finds itself at an infinitely remote 
point. If the values > l/Sap are possible in the PGW metric ([Ml) , it 

makes no sense to continue the solution of magnetohydrodynamics beyond the 
hypersurface u = u^,, since at the shock wave front all the invariants suffer a 
secondtype discontinuity. 

6 Source of energy 

Since on a singular hypersurface u — u^, the energy densities of the plasma and 
the magnetic field tend to infinity, while the velocity of plasma as a whole tends 
to the speed of light, the total energy of a MHD shock wave tends to infinity. 
Thus, as a result of the drift, the the complete plasma energy grows (indeed, 
it grows to infinity), and it is necessary to reveal the nature of a source of this 
energy. Note that the MHD equations were solved against the background of 
the PGW metric. In this sense the PGW is an inexhaustible source of energy 
for the magnetoactive plasma; it is precisely for this reason that a singular state 
is formed. A singular state originating in the plasma under the infiuence of 
PGW, disturbs the basic assumption ([T]) on the weakness of the GW-plasma 
interaction. In a more comprehensive self-consistent problem with regard to 
gravitation the allowance made for the impact of the shock wave on the PGW 
should lead to a PGW amplitude damping to the values 



Thus the discovered effect of MHD wave generation by a gravitational wave 
(GMSW) can be an effective mechanism of GW energy absorption. The author 
is unaware of any other GW energy absorption mechanism of similar efficiency. 

A strict solution to the problem of PGW energy transformation to the energy 
of a shock wave is possible only on the basis of a study of the self-consistent 
set of the Einstein and MHD equations. However, a qualitative analysis of the 
situation may be performed using a simpler model of the energy balance. Let 
us consider the case of strictly transverse PGW propagation {ft = tt/2) in a 
nonrelativisic plasma. Then the energy fiux of magnetoactive plasma should 
have the direction of PGW propagation, i.e., along the Ox^ axis. Let f3o{u) be 
the PGW vacuum amplitude and /3(m) be its amplitude with allowance made 
for absorption in the plasma. From summed energy fiux conservation it follows: 



max(l^l) < l/2a^. 



(181) 



T^\(3)+t^\l3)^t^\(3o), 



(182) 
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where t^^ is the energy flux of a weak GW in the direction Ox^ (see [11]). Using 
the solution to the MHD equations for nonrelativistic plasma (|174p in the case of 
strictly transverse PGW propagation, as well as the expressions for the summed 
plasma EMT and weak PGW energy flux [11], we reduce (I182p to the form 
(returning to the standard system of units) 



ttGeq 



where /3 means the derivative with respect to the time t. Setting in what follows 

(w being the GW frequency), ctg ^ 1 and ^ 1 and tending A to zero, we 
reduce the latter equation to the form: 



X 



(l-2ag/3)4 
where x is a dimensionless parameter: 



(183) 



(184) 



^TiGE^jc?. The approximation ([T]) is equivalent to the condition: 

« 1- 



(185) 



In the conditions of strong PGW energy absorption we have ^ /Sg, therefore 
in this case (|184|) implies: 



1/2 



/3 = 



2al 



X 

l/3o| 



1/2- 



(186) 



It follows from (|186p that the conditions for a sufHciently strong absorption of 
a weak PGW are: 



> 1 ; 



1 



3o| >X. 



(187) 



The PGW amplitude damping factor 7 is: 

\0\ 1 



7 



Wo\ 2al\l3o\ 



1 - 



X 



1/2 



(188) 



7 Astrophysical consequences of GMSW 

There naturally arises a question, whether or not the GMSW exist in nature, i.e., 
can the condition for their creation (|149l) be realized? In laboratory experiments 
the above necessary condition (2q:^|/3| > 1) is unattainable due to the smallness 
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of both the gravitational waves in terrestrial conditions and stable laboratory 
magnetic fields. For galactic magnetic fields {H ~ fO~^ fO^^ Gauss) and 
interstellar medium (p ^ lO^^^g/cm^) a'^ ^ 10^* ), so GMSW in the 
interstellar medium do not arise. 

The only possible source of a GMSW can apparently be neutron star magneto- 
spheres on the stage of quadrupole neutron oscillations, as well as on stage 
of a Supernova. Table 1 represents the results of calculating the parameter 
2a^(|/3|), ((1/3 1) is the average amplitude of irradiated GW) for neutron stars 
magnetospheres, performing quadrupole oscillations of the fundamental quad- 
rupole mode (n = 0). The data placed in the first four columns are taken (or 
recalculated) from the book [12]. The average GW amplitude is calculated us- 
ing standard formulas, from the average gravitational radiation power, N, (see 
for instance [11]) under the assumption that the average relative amplitude of 
quadrupole oscillations of a neutron star (the parameter ((5i?/i?)^)^/^ in the 
book [12]) is equal to 10^®. It is also assumed that the neutron star magneto- 
sphere consists of completely ionized hydrogen. In addition, the parameter uq 
(electron number density on the neutron star surface) was found from the con- 
dition T = 1 , where t = J nadl is the optical thickness of the magnetosphere. 
It is assumed also, that the scale of dense magnetosphere is of the order of the 
neutron star radius R. 

Note that the electron concentration values near a neutron star surface, 
obtained from the condition that the optical thickness of the magnetosphere 
equals unity, is approximately by 4 orders of magnitude larger than the estimates 
of no from [19], obtained on the basic of a dimensional analysis of the Maxwell 
equations. If we accept Pacinis estimates for no [19], then the GMSW factor 
will grow by 4 orders as compared with that of Tabic 1 . 

Table 10 The GMSW parameter: 2a^{\l3\) in a neutron star magnetosphere 



M/Mq 


R 




N/m/RY) 


m) 


no 


2a-{m 


0,405 


5,00 


5249 


1,2(50) 


2,93(-10) 


5,97(18) 


0,260 


0,682 


8,42 


2,02(4) 


2,9(53) 


4,20(-9) 


3,54(18) 


3,32 


0,677 


12,60 


8987 


7,0(52) 


l,64(-9) 


2,36(18) 


3,670 


1,954 


9,99 


1,66(4) 


1,6(55) 


l,41(-8) 


2,99(18) 


30,0 



M/Mq is the neutron star mass related to the Solar mass; R is the star radius in 
km.; oj is the gravitational radiation frequency in sec" N/{{5R/RY) is the stars 
gravitational radiation power in erg/sec; {|/3|) is the GW average amplitude; no is the 
concentation of electrons on the stellar surface in cm~^ . 



For a star with a mass of 0,682 Mq we get from Table 1: % = 1,33 • 
10""; X/{\M) = 0,317. According to ([TMl) ) we find the GW damping fac- 
tor: 7 = 0, 206. Thus, the GW amplitude in this case is damped by a factor 

^quantities in paranthcses mean the degrees of 10 
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of 5 and the gravitational radiation power by a faetor of 25. For a star with a 
mass of 1, 954Mq x = 1, 42 • 10""; x/(|/3o|) = 1,01- IQ-^ -^ 7 = 0, 032, i.e., in 
this case the GW amphtude can be diminished by a factor of 30 times and the 
GW power by a factor of 900! So in neutron star magnetospheres there exist 
necessary conditions for GMSW excitation. 

If the magnetic field of a neutron star is described as that of a dipole, then the 
geographic angle (counted relative to the magnetic equator) will be connected 
with the above angle by the relation = 7r/2 — 0. Therefore the GMSW 
excitation condition depends on the angle 9: 

sin^ e < 1 - ^ 



Thus, in the magnetosphere of a neutron star (or a Supernova) a GMSW 
can be excited in the region of the magnetic equator, similarly to pulsars with 
a knife-like radiation pattern. In this region, as was demonstrated by the above 
examples, the gravitational radiation can be absorbed practically completely 
by the excitation of shock waves. A neutron star of such type should radiate 
gravitational waves only from its magnetic poles, similarly to pulsars with a 
pencil-like radiation pattern. In this case the probability of observing a GW 
sources can be sharply dropped. However, suddenly the GMSW open another 
way of observing gravitational waves. If such a shock wave is formed and takes 
off the magnetosphere, it will carry with itself (at a subluminal speed) super- 
strong magnetic fields into the interstellar space. The interaction of cosmic 
plasma with such magnetic fields can lead to the anomalous electromagnetic 
phenomena in the radio and optical spectral ranges. It should be stressed that 
there is no other mechanism able to accelerate a shock wave to subluminal 
velocities. 

Thus, the GMSW phenomenon may displace the centroid of GW experi- 
ments from direct GW detection to optical observations of the Supernovae and 
their remainders. However, to be quite certain that the GMSW do exist in 
neutron star magnetospheres and to know the observational manifestations of 
this phenomenon, it is necessary to solve a number of problems: 

1. To solve the self-consistent problem of GW propagation in a homogeneous 
plasma with an embedded magnetic field and calculate the damping dec- 
rement; 

2. To clarify the effects on the this phenomenon from plasma inhomogeneity 
and external gravitational fields; 

3. To clarify the effect on this phenomenon from the curvature of magnetic 
lines of force; 

4. To investigate optical manifestations of GMSW. 
We intend to study these problems in future papers. 
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